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[i]  Differential  cross  sections  for  Raman  scattering  between  the  multiplet  levels  of  the 
ground  electronic  state  of  atomic  oxygen  are  calculated  from  first  principles.  The 
calculated  cross  sections  are  in  order-of-magnitude  agreement  with  the  available 
experimental  measurements.  The  use  of  a  Raman  LIDAR  (Light  Detection  and  Ranging) 
has  the  potential  of  giving  a  three-dimensional  map  of  the  density  of  each  of  the  three  fine 
structure  levels  of  atomic  oxygen.  Since  it  was  established  earlier  that  the  fine  structure 
levels  are  in  local  thermodynamic  equilibrium  (LTE)  up  to  at  least  350  km  altitude,  the 
Raman  LIDAR  therefore  also  has  the  potential  of  giving  a  three-dimensional  map  of  the 
temperature  of  the  atmosphere.  INDEX  TERMS:  0310  Atmospheric  Composition  and  Structure: 

Airglow  and  aurora;  0355  Atmospheric  Composition  and  Structure:  Thermosphere — composition  and 
chemistry;  0340  Atmospheric  Composition  and  Stnicture:  Middle  atmosphere— composition  and  chemistry; 
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oxygen,  Raman  transitions  between  fine  stnicture  levels,  Raman  LIDAR,  local  thermodynamic  equilibrium, 
three-dimensional  map  of  atmospheric  temperature  and  atomic  oxygen  density 
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1.  Introduction 

[2]  The  role  of  atomic  oxygen  in  the  heat  budget  of  the 
terrestrial  thermosphere  can  hardly  be  overestimated 
[Sharma  and  Rohle ,  2002],  Three  important  processes 
involving  atomic  oxygen,  namely  the  15  pm  emission 
pumped  by  the  impacts  of  O  with  C02,  5.3  pm  emission 
pumped  by  the  impacts  of  O  with  NO,  and  63  pm  fine- 
structure  emission  O(3P0  — ►  0(3P2)  are  the  dominant 
mechanisms  for  cooling  the  thermosphere  and  determin¬ 
ing  its  temperature  and  density  structure  [Sharma  and 
Rohle ,  2001].  It  is  therefore  very  desirable  to  know  the 
temperature  and  density  of  atomic  oxygen  in  the  thermo¬ 
sphere.  Recent  satellite-based  passive  remote  sensing  of 
atomic  oxygen  density  in  the  130-175  km  altitude  range 
by  measuring  the  63  pm  limb  emission  by  Cryogenic 
Infrared  Spectrometers  and  Telescopes  for  the  Atmosphere 
(CRISTA)  aboard  ASTRO-SPAS  (Astronomical  Shuttle 
Pallet  Satellite)  for  the  November  1994  to  August  1997 
period  [Grossmann  et  al ,  2000]  have  given  very  surpris¬ 
ing  results:  “The  resultant  densities  are  on  the  average 
about  40%  below  the  predictions  of  the  MSIS  model  for 
the  low  solar  activity  conditions  of  the  CRISTA  experi¬ 
ments  and  at  all  altitudes  covered  by  the  measurement. 
The  densities  exhibit  considerable  horizontal  structures 
which  are  only  in  partial  agreement  with  the  MSIS 
predicted  variations.  The  latitudinal  distribution  shows 
steeper  decreases  towards  the  poles  than  given  in  MSIS.” 

[3]  The  CRISTA  measurements,  while  showing  consid¬ 
erable  horizontal  structure,  also,  because  they  are  limb 
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measurements,  average  the  structure  over  the  line-of-sight 
and  partially  mask  it.  Since  it  is  possible  to  gate  the 
scattered  Raman  pulse,  one  can  determine  the  O  density 
and  temperature  in  a  particular  space-time  interval,  giving 
us  a  method  of  determining  these  two  very  important 
parameters  for  atmospheric  modeling  without  limb  averag¬ 
ing  over  a  couple  of  thousand  kilometers.  It  is  hoped 
that  this  work  will  be  as  useful  to  future  experiments  as 
our  previous  work  [Sharma  et  al ,  1988,  1990,  1994; 
Zachor  and  Sharma ,  1989]  was  to  CRISTA  experiments. 
Assuming  the  fine  structure  levels  of  atomic  oxygen  to  be  in 
local  thermodynamic  equilibrium  (LTE),  Sharma  et  al. 
[1988]  showed  that  the  63  pm  and  145  pm  transitions 
could  be  inverted  in  the  90-250  km  altitude  range.  The 
simulated  signal  with  noise  incorporated  was  inverted  in 
the  145  pm  channel  by  Zachor  and  Sharma  [1989]  and 
in  the  63  pm  channel  by  Sharma  et  al.  [1990].  Since  the 
idea  that  the  fine  structure  levels  were  in  LTE  was  not 
generally  accepted,  the  newly  calculated  cross  sections  for 
transitions  between  the  fine  stnicture  levels  during  0-0 
impacts  [Zygelman  et  al. ,  1994a,  1994b]  were  used  to 
demonstrate  that  these  levels  were  always  in  LTE  at  least 
up  to  350  km  altitude  [Sharma  et  al ,  1994]. 

[4]  Atoms  with  half-filled  or  completely  filled  outer 
electronic  orbit  have  only  isotropic  polarizability  and 
therefore  do  not  Raman  shift  the  frequency  of  the 
incident  radiation.  Atoms  with  neither  half-filled  nor 
completely  filled  outer  electronic  orbit,  on  the  other  hand, 
possess  both  isotropic  and  anisotropic  polarizabilities  and 
therefore  do  Raman  shift  the  frequency  of  the  incident 
radiation.  From  these  considerations,  Schlossberg  [1976] 
pointed  out  that  the  fluorine  atoms  might  be  detected  by 
monitoring  the  Raman  scattered  radiation  and  estimated 
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that  about  1016  atoms  could  be  detected  by  this  method. 
The  atomic  oxygen  ground  electronic  state  (3P)  has  eight 
electrons  ls22s22p4  and  its  outer  electronic  orbit  is  neither 
half-filled  nor  completely  filled.  It  therefore  has  nonzero 
anisotropy  of  polarizability;  hence  Raman  scatters  inci¬ 
dent  radiation.  Dasch  and  Bechtel  [1981]  measured,  with 
an  excitation  source  of  532.1  nm,  the  Raman  scattering 
cross  section  without  change  in  polarization  of  ground 
state  atomic  oxygen.  These  authors  also  calculated  the 
Raman  scattering  cross  section  assuming  that  the  3Si 
electronic  state  located  at  76,700  cm'  1  is  the  only  state 
that  makes  a  significant  contribution  to  the  Raman  cross 
section.  Since  Dasch  and  Bechtel  [1981]  do  not  give 
details  of  their  calculation,  and  since  the  transition  mo¬ 
ment,  derived  by  us  in  Appendix  A,  differs  by  a  factor  of 
1.4  from  that  given  by  these  authors,  we  repeat  their 
calculation  using  the  full  power  of  the  angular  momentum 
algebra  [ Zare ,  1987].  We  obtain  results  that  are  in  order- 
of-magnitude  agreement  with  the  available  measurements. 

[5]  Section  2  gives  the  formulation  of  the  equations;  the 
results  are  computed  and  compared  with  the  earlier  mea¬ 
surement  and  calculation  of  Dasch  and  Bechtel  in  section  3. 
Section  4  states  the  possible  application  of  these  results  to 
measuring  the  atomic  oxygen  density  and  temperature  in  the 
terrestrial  thermosphere. 


unit  frequency  interval  per  unit  volume  in  solid  angle  dD  is 
given  by 


M£)  =  ^ Ldo. 

(2tvc)‘ 


(4) 


It  is  usual  in  the  literature  to  incorporate  an  additional  factor 
of  two  in  equation  (4),  corresponding  to  the  two 
independent  polarization  directions  of  the  emitted  photon. 
Since  we  calculate  the  scattering  cross  section  for  the  two 
directions  of  polarization  separately,  our  density  of  states 
will  not  have  that  additional  factor  of  two. 

[7]  We  make  the  dipole  approximation  and  assume  that 
the  dimensions  of  the  atoms  and  molecules  we  are  dealing 
with  are  veiy  short  compared  to  the  wavelength  of  the 
electromagnetic  radiation  and  approximate  exp(ik.r)  term  in 
equation  (3)  by  unity.  Using  the  relationship  between  the 
expectation  values  of  the  momentum  and  dipole  operators 

-(n|p|m)  =  kx)mil<n|p|m),  (5) 

F 

where 


2.  Formulation 

[5]  The  transition  probability  for  going  from  state  a  to 
state  b  with  the  absorption  of  a  photon  of  frequency  10  from 
the  electromagnetic  field  and  the  emission  of  one  of 
frequency 

{Eh  ~  Ea) 

Wba  =  - - 7 - 


in  solid  angle  dft  is  given  by  Fermi’s  Golden  Rule 
[ Kramers ,  1964] 

P  =  ^  /  dio|<b|W(w,t)|a)|2&(u)ba  -  w)p(w)da  (1) 

where  p(w)  is  the  density  of  states  per  unit  frequency 
interval.  W(oj,  t)  is  the  interaction  energy  of  the  electro¬ 
magnetic  field  of  frequency  uj,  and  a  spinless  particle  of 
mass  p,  momentum  p,  and  charge  e  given  by 

W(w,t)  =  — —  (A.p)+-%A2  (2) 

pc  2pc- 


A(uj,  t)  =  e(exp(— i(k.r-wt)  +exp(i(k.r  -  ut)).  (3) 


The  electromagnetic  field  is  normalized  to  contain  one 
photon  of  frequency  oj  in  unit  volume  and  e  is  a  unit 
(polarization)  vector  in  the  direction  of  the  electric  field 
perpendicular  to  the  propagation  vector  k  of  the  radiation. 
The  density  of  states  p(oo)  of  the  electromagnetic  field  per 


is  the  transition  frequency  between  levels  n  and  m.  We  can 
now  write  the  first-order  matrix  elements  of  the  interaction 
operator  W  between  states  n  and  m 

(n|Wjm)  =  —  iw„m  (^  (eiu*  +  e"“l) .  (6) 

Raleigh  and  Raman  scattering  are  second-order  processes, 
i.e.,  in  going  from  the  initial  state  a  to  final  state  b  the 
interaction  operator  acts  twice.  The  second-order  interaction 
operator  consists  of  three  tenns:  (1)  An  incident  photon  of 
frequency  u;  is  absorbed  and  the  particle  (atom  or  molecule) 
undergoes  a  transition  to  virtual  state  |v)  followed  by  the 
emission  of  a  photon  of  frequency  ujs  with  the  particle  going 
to  the  final  state  |f ).  (2)  A  photon  of  frequency  ios  is  emitted 
and  the  particle  undergoes  a  transition  to  virtual  state  |v) 
followed  by  absorption  of  a  photon  of  frequency  lj  with  the 
particle  going  to  the  final  state  |f).  (3)  A  photon  of 
frequency  io  is  simultaneously  emitted  and  absorbed  and  the 
particle  stays  in  the  initial  state. 

[s]  The  first  two  tenns  are  called  the  Kramers-Heisenberg 
terms  while  the  third  one  is  called  the  Waller  term 
[. Mizushima ,  1970].  The  Waller  term  is  also  called  the 
A2  tenn  because  it  arises  due  to  this  interaction  in  the 
perturbation  potential  (equation  (2)).  An  interesting  point 
deserves  mention  here:  The  atom  acts  purely  as  a  spectator 
in  the  Waller  or  the  A2  term,  i.e.,  undergoing  no  change  at 
all.  This  tenn,  although  4-5  orders  of  magnitude  smaller 
[Loudon,  1983],  is  still  needed  to  obtain  the  commonly 
given  expression  for  Rayleigh  and  Raman  scattering  given 
below  (equation  (7)). 
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[9]  The  three  terms  summed  together  give  the  differential 
cross  section  for  scattering  photons  [Mizushima,  1970; 
Kramers ,  1964;  Sharma  and  Levin,  1973]: 

da  /(fl(es.n,)|v)(vie.|i|i)  (fl(e.M.)|v)(v|es.Hi)>\  ' 

dQ_  c4  VI  (Ev-Ej-M  (Ev-Er  +  M  } 

(7) 


where 

hu  is  the  energy  of  the  incident  photon 
fiu;s  is  the  energy  of  the  scattered  photon 
|i)  is  the  wavefunction  of  the  initial  state  of  the  particle 

|v)  is  the  wavefunction  of  the  virtual  state  of  the  particle 

jf)  is  the  wavefunction  of  the  final  state  of  the  particle 
\i  is  the  dipole  moment  operator 
e  is  the  polarization  (electric)  vector  of  the  incident 
photon 

es.  is  the  polarization  (electric)  vector  of  the  scattered 
photon. 


[11]  The  ground  state  of  atomic  oxygen  is  3Pj,  j  =  1,2, 
and  3.  The  different  j  levels  are  split  because  of  the  spin- 
orbit  interaction  with  3P2  being  the  lowest  at  0  cm"1, 3Pi  at 
158.2687  cm"1,  and  3P0  at  226.9851  cm"1  [Zink  et  al , 
1991].  Using  the  Wigner-Eckart  theorem,  the  matrix  ele¬ 
ment  of  e.\i  between  the  initial  and  virtual  states  can  be 
written  [. Zare ,  1987] 


/A  \  1  „ 

(n" £" s"j " m" | e . (i | nf sj m )  =  ( yj‘(n'r^iK)(-l)s+1_M  K" 

x  [(2r  +  l)(2j+l)pW(^"j";sl) 
x  (('"  ||  Y,  ||  £)C(j<j";m,0,m)6  mm" 

(ii) 


and  the  matrix  element  in  the  uncoupled  representation 


[t 


Y, 


||f)  =  C(Gf";000) 


'  3(2f  +  1)  ‘ 
4tv(2  f  +  1)_ 


(12) 


[10]  We  now  write  [Zare,  1987] 

1 

£'1JL=I  (8) 
M=—  I 

where  p0  (£0),  E  i  (£-i)»  Ei  (£1)  are  projections  of  the 
vector  |x  (e)  on  the  unit  vectors  along  the  z,  (l/2)1/2(x  —  iy), 
and  — (l/2)1/2(x  +  iy)  axes,  respectively.  The  projection 
8,(e_i)  corresponds  to  left  (right)  circular  polarization. 
Radiation  polarized  along  x-axis  is  obtained  by  the  coherent 
(summing  before  taking  absolute  square)  combination 
(l/2)1/2(e_i  -  £1)  while  that  polarized  along  y-axis  is 
represented  coherently  by  the  combination  (i/(2)1  “)(£_]  + 
z{).  Recalling  that  pm  =  (4ir/3)1/2  pYlni(Qft)  we  can  write 


giving 

<nVV7m'>.p|nfsjm>  =  <n"f"|p|n e)C{£U";  00)  (-)s+1  bss» 

x[(M+l)(2j+l)]*W(fj<"j";s0 

x  C(j lj;/; m, 0, m)8mmw  (13), 

where  {n"£"|p|nf)  is  the  radial  matrix  element  of  the  dipole 
moment  p,  C  are  the  Clebsch-Gordan  coefficients,  and  W  is 
the  Racah  coefficient.  The  third  “magnetic  or  M”  quantum 
number  in  a  Clebsch-Gordan  coefficient  is  equal  to  the  sum 
of  the  first  two  and  is  omitted  when  there  is  no  danger  of 
confusion.  The  matrix  element  of  es.[i  between  the  virtual 
and  final  states  can  be  similarly  written  as 


(4  M=+l 

“Tj  P  (9) 

5  /  M=-l 

Similarly,  we  can  write 

m=+l 

T U  Y,  (io) 

To  derive  the  scattering  angle  and  polarization  dependence 
of  the  emitted  radiation,  we  follow  Sharma  and  Levin  [1973] 
and  consider  two  coordinate  systems:  one  fixed  in  the 
incident  beam  (coordinate  system  1 )  and  the  other  fixed  in 
the  scattered  beam  (coordinate  system  2).  We  take  the  y-axes 
of  the  two  beams  along  their  respective  directions  of 
propagation.  The  z-axis  of  the  two  coordinate  systems  is 
taken  perpendicular  to  the  plane  containing  the  incident  and 
scattered  beams.  The  two  coordinate  systems  thus  have  the 
same  z-axis  and,  in  the  experimental  setup  of  Dasch  and 
Bechtel  [1981],  coordinate  system  1  can  be  transformed  into 
coordinate  system  2  by  a  90°  counter  clockwise  rotation 
about  this  axis.  Further,  the  incident  beam  is  polarized  along 
the  z-axis.  Thus  we  are  restricted  to  M  in  Equation  9  to  only 
one  value,  M  =  0.  Further,  without  loss  of  generality,  we  can 
pick  our  space-fixed  coordinate  system  to  be  coincident  with 
the  coordinate  system  2. 


(nY s'j'm'|es.jx|n'Y,s"j"m")  =  tfi'W't)  £  (-)%£q 

q 

x  C( j'Vj';  m",  — q,  m"  —  q) 
x  (— 0s"+1_r_r&S's"[(2/'  +  1) 
x  (2j"+  l)]-W(^"j'Yj';  s'V) 
x  C(C"lC';00)  (14), 

where  eq  is  a  unit  vector  along  the  spherical  axis  and  cq  is 
the  projection  of  the  unit  vector  along  the  polarization 
direction  of  the  scattered  beam  in  the  coordinate  system 
number  2,  the  space-fixed  coordinate  system. 

[12]  Using  the  relations 

C(jlj";  m.  0,  m)C(j"lj';  m,  -q,  m  -  q) 

=  Y  I(2j"  +  l)(2f  +  l)]‘W(jlj'l;j"f) 

f 

x  C(1  If; 0,  — q,  — q)C(jfj';  m,  — q,m  —  q)  (15) 

and 

C(jfj';m,-q.m-q)  =  C(jj'f;m, -m  +  q, q) 

(16) 
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we  can  write 

(f|(es.|i.)|v)(v|ej.|i|i)  =  {n'f|ii|nY'){n"^"||i|n/,)6ss«6s<S” 

x  m[(2j+  l)(2j"+  l)2 

x  (2j'+  l)(2f"+  l)(2f  +  l)]!C(fir;000) 

x  C{("  I C:  000)W((jrj";sl)W(/"j"('j';  s"  i ) 

x  H  (— )qcqeq  ]T  c( jj'f;  m,  -m  +  q.  q) 
q  f 

x  C(1  If;  0,  — q,  -q)W(jlj,l;j,/f)  (17) 


If  the  scattered  radiation  has  the  same  polarization  as  the 
incident  beam,  i.e.,  it  is  polarized  in  the  z  direction,  then 
q  =  o,  co  =  1,  and  f  =  0,  2.  Since  f  has  to  be  nonzero  for 
transition  from  one  multiplet  level  to  another,  only  f  =  2  is 
allowed.  The  combination  j  =  1,  j'  =  0,  and  f  =  2  does  not 
satisfy  the  triangular  inequality  leading  to  the  Racah 
coefficient  W(1 101  ;1 2)  =  0;  thus  the  transition  3P,  -»  3P0, 
shifted  68.7165  cirT1  from  the  incident  frequency,  without 
change  in  the  direction  of  the  polarization  vector  is  not 
Raman  allowed.  The  differential  cross  section  for  Raman 
scattering  from  multiplet  level  j  to  level )  without  change  in 
the  polarization  is  then 


Averaging  over  the  initial  magnetic  quantum  numbers  m 
and  summing  over  the  final  ones  for  a  given  value  of  q  we 
get 

^<j  - n  =  ^  £  m  +  ')C!("f:0. -q. -q)|c,l' 

S  ((E,  -  Ei  -  M  +  (E,  -  E,  +  m) 

x(n7'||i|n"f")  x  (n"/’"||i|nf> 

x  x(_,)2.-f-w*-/'|-(2j"+  i)2(2r  +  l)(2f  +  1)]I6S'S6S»'S 
xC(f  1  f";  000)C(f"  if ;  000) 
x  W(fjf"j";  si  )W(f"jYj';  s"  1  )W(  j  I  j'  1  :j"f) 

(18) 

In  equation  (18)  f  denotes  the  rank  of  the  polarizability 
tensor  and  q  denotes  the  direction  of  the  polarization  vector 
of  the  scattered  photon. 

3.  Computation  of  Raman  Scattering 
Cross  Sections 

[13]  Equation  (18)  is  used  to  calculate  Raman  Scattering 
cross  sections  between  the  multiplet  components  of  atomic 
oxygen  with  different  j  values.  The  Kronecker  deltas  require 
that  the  virtual  states  must  also  be  spin  triplets.  Also,  £=£'=  s  = 
s'=  1.  For  the  Racah  coefficient  W(  jlj'l;j"f)  with  j  ^  j'to  be 
nonzero,  the  index  f  must  be  nonzero.  Therefore  only  aniso¬ 
tropic  components  of  the  polarizability  tensor  can  Raman 
scatter  between  the  multiplets  of  the  ground  state  of  atomic 
oxygen .  Following  Dasch  and  Bechtel  [  1 98 1  ],  we  assume  that 
only  one  virtual  state,  the  3S)  state  (76700  cm-1),  makes  a 
significant  contribution  to  the  scattering  process.  Noting  that 
t  =  0  and ]n  =  1  gives  W(lj01;ll)  =  W(011j/;11)=  1/3.  We 
also  note  that  £  =  t  =  1,  giving  C(110;00)  =  — ( 1/3)1/2  and 
C(01 1  ;00)  =  1 .  Substituting  these  values  as  well  as  the  dipole 
moments  from  Appendix  A,  we  get 


-/)  =  ^  (L2  X  (2j'  +  1)C2(1 12;  0, 0,0) 


1 


^(Ev  -  Ei  -  M  (Ev  -  Ef  +  M 
x  W(jljM;  12) 


(20) 


For  the  incident  radiation  of  532.1  nm  wavelength,  the 
differential  cross  section  for  the  3P0  ->  3P2  Raman  transition 
is  2.0  x  10-31  cm2/str  and  that  for  the  3Pi  ->  3P2  Raman 
transition  is  1.5  x  10“31  cm2/str;  the  ratio  of  the  two  cross 
sections  is  just  the  ratio  of  the  squares  of  the  Racah 
coefficients 


/W(  1121;  12)\  2_  3 
\ (0121;  12) /  ~  4 

The  experimentally  measured  cross  sections  are  (4  ±  1)  x 
10-31  and  (6  ±  1)  x  10-31  cm2/str,  respectively.  The 
calculation  of  Dasch  and  Bechtel  [1981]  gives  1.2  x 
10-31  cm2/str  for  3P0  ->  3P2  and  2.8  x  10-ri  cm2/str  for 
3Pj  3P2.  The  ratio  R  of  the  cross  sections  calculated  by  us 
3P{  3P2/3P0  — y  3P2  =  0.75  while  that  obtained  by  Dash 

and  Bechtel  is  2.3;  the  experimental  ratio  being  1.5  ±  0.6. 

[14]  The  expression  for  the  differential  cross  section  for 
Raman  scattering  with  change  in  the  direction  of  polariza¬ 
tion  is 


jg(J— f)=  £  «'+» 

C  y  f=1.2  q  — —  l.+  I 

x  C2(l If; 0,  — q,  — q)|cq|2 
(  1  1 
X  \(EV  -  Ej  -  i 


■  + 


hut)  (Ev  —  Ef  +  fiu;) 


x  W(jlj'l;lf) 


(21) 


For  the  x  direction  of  the  polarization  vector  of  the  scattered 
photon,  q  =  ±1  and  |c_,  |2  =  |c,  |2  =  1/2.  Using  C2(l  1 1  ;0  -  1 )  = 
C2(l  1 1  ;01)  =  C2(l  12;0,- 1)  =  C2(l  12;01 )  =  1/2,  we  can  write 
the  expression  for  the  differential  cross  section  as 


da 

dQ 


(j 


uw;1  (1.2  x  10 -,8)4 

J  /  —  n 


da 


^(2j'  +  l)C2(llf;0,-q,-q)  —  (j 


r.q 


X  C„  X 


-  +  - 


^Ev  -  Ej  -  h w)  (Ev  -  Er  +  M 
xW(jlj'l;lf) 


(19) 


=  l0~lg)-lx  V(2j'+1) 

c  9  2  f"2  v 

/I  1 

X  V(Ev  -  Ei  -  *w)  +  (Ev  -  Ef  +  fiu;) 

x  W(j]j'l;  lf)| 


(22) 
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Table  la.  Differential  Cross  Section  (cm2/str)  Raman  Scattering 
for  532.1  nm  Radiation  (xlCT31)  Without  Change  of  Polarization 


Transition 

Frequency 
Shift,  cm-1 

Present 

Calculation 

Dasch  and 
Bechtel  [1981] 
Calculation 

Dasch  and 
Bechtel  [1981] 
Experiment 

3Po  -  3P2 
3Pt  3P-> 
3Po  -  3Pi 

226.9851 

158.2686 

68.7165 

2.0 

1.5 

0.0 

1.2 

2.8 

0.0 

4±  1 

6  ±  1 
<1 

where  we  have  written  ooba  =  eo  and  e  is  the  polarization 
vector  of  the  emitted  photon  and  \i  is  the  dipole  moment 
operator.  Anticipating  that  the  absolute  square  of  the  matrix 
element  will  be  shown  to  be  angle  independent,  we  integrate 
over  the  solid  angle  obtaining 

p=T"1=?03|e'(bWa)|2-  (A3) 


The  f  =  1  term  gives  zero  for  the  differential  cross  section  of 
the  3P<)  — *  3P2  transition,  since  W (0 1 2 1 ;  1 1 )  =  0.  The 
differential  cross  section  of  the  3P0  ->  3Pi  Raman  transition, 
which  is  zero  when  the  polarization  is  unchanged,  is  1.5  x 
KT31  cm2/str.  The  cross  section  for  the  3Pj  — >  3P2  transition 
is  0.6  x  10“31  cm2/str. 

[15]  The  f  =  2  term  gives  zero  for  the  differential  cross 
section  of  the  3P0  — ►  3Pi  transition  because  W(0 1 1 1 ;  12)  =  0. 
The  cross  sections  with  z  to  x  change  of  polarization  for  the 
3Pj  — >  3P2  and  3P0  — >  3P2  transitions  are  1.1  x  10“31  and 
1.5  x  10~31  cm2/str,  respectively.  The  differential  cross 
sections  for  Raman  scattering  with  z  to  x  change  of 
polarization  are  1.5  x  10-31,  1.5  x  10  31,  and  1.7  x 
10~31  cm2/str  for  the  3Po  — 3 -  3P2>  3Po  — »  3Pi,  and  3Pi  — >  3P2 
transitions,  respectively.  These  results  are  summarized  in 
Tables  la  and  lb. 

4.  Discussion 

[16]  We  have  calculated  Raman  scattering  cross  sections 
among  the  multiplet  fine  structure  levels  of  atomic  oxygen. 
We  have  shown  that  the  isotropic  component  (tensor  of  rank 
zero)  of  the  polarizability  tensor  makes  no  contribution  to 
this  cross  section  and  have  calculated  the  contributions  of 
the  polarizability  tensors  of  rank  one  (vector)  and  two.  This 
calculation,  as  well  as  that  of  Dasch  and  Bechtel  [1981], 
gives  order-of-magnitude  agreement  with  the  available 
experimental  data.  Considering  the  approximation  of  taking 
into  account  only  one  virtual  state,  this  agreement  cannot  be 
considered  unreasonable.  The  Raman  LIDAR  experiment 
proposed  here  would  permit  a  three-dimensional  map  of 
each  of  the  three  fine  structure  components  of  the  ground 
electronic  state  of  atomic  oxygen.  Because  these  fine 
structure  levels  are  in  LTE  up  to  at  least  350  km  altitude 
[Shanna  et  al ,  1994],  a  three-dimensional  map  of  these 
levels  would  also  give  a  three-dimensional  map  of  the 
temperature. 

Appendix  A 

[17]  To  evaluate  the  radial  integral  from  the  radiative 
lifetime  of  the  3Sj  3Pj,  J  =  0,1,2,  we  start  with  Fermi’s 
Golden  rule 

P  =  T-'  =|£  J dulOjIWKOlaJl^tuta-uJpMdn,  (Al) 

where  t  =  4.8  ns  [ Wiese  et  al. ,  1966]  is  the  radiative  lifetime 
of  the  3Sj  — >  3Pj  transitions  at  76,700  cm"1,  and  the  matrix 
element  is  the  transion  dipole  moment  between  the  upper 
and  lower  level.  Proceeding  as  in  section  2,  we  can  write 

P  =  T~'  =  2t^  (c)  /  dn'S  '  ^bWa^2’  (A2) 


Without  loss  of  generality  we  take  our  space-fixed  z-axis 
along  the  direction  of  propagation  of  the  emitted  photon. 
Using  the  Wigner-Eckart  theorem,  we  obtain 

e  •  (n'fs'j'mVKsjm)  =  +  l)(2j  +  1)] 

x  W(fj^j';sl)(nY'  ||  Y,  ||  nl)<nY|nK> 
x  Y  e-MC(jlj';m,M,m')  (A4). 

Using  the  relations 


(n't  ||  Y,  ||  n£)  =  C(e\t;00) 


'  3(2 i+  1)  ' 
4ir(2/,'  +  1). 


and 


C(jlj';m.M.m  +  M)  =  (-iy_m 


x  C(jj'l;  m,  -m  -  M,  -M), 


(A5) 


(A6) 


we  get 

e  •  (nYs'j'mVKsjm)  =  C(Of;000)(nT|Hn^(-rH'+j~m6Ss' 

..  r(2^+l)(2j+l)(2j'  +  l)li 

[  3 

xW(^j';sl)  Yj  £-m 

M=-l.+  l 

x  C  ( j  j '  1 ;  m,  —  m  -  M,  — M).  (A7) 

Taking  the  absolute  square,  averaging  over  the  initial 
magnetic  quantum  number,  and  summing  over  the  final 
ones  and  over  M  also,  we  get 


|£-(fs'j'|pKsj)|2-C2(flf;000) 


2(21  +  offi  +  iy 

3 


x  W2(fjfj';  sl)|{nY|p|n£)|2 


(A8) 


Table  lb.  Differential  Cross  Section  (cm2/str)  Raman  Scattering 
for  532.1  nm  Radiation  (xl(T31)  With  Change  of  Polarization 


Transition 

Frequency 
Shift,  cm-1 

Present 

Calculation 

Dasch  and 
Bechtel 
[1981] 
Calculation 

Experiment 

3p(,  —  3p, 
3p,  —  3p. 
3Po  -  ’Pi 

226.9851 

158.2686 

68.7165 

1.5 

1.7 

1.5 

0.9 

1.4 

0.0 

Not  available 
Not  available 
Not  available 
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Recalling  that  the  initial  state  is  3Sj,  i.e.,  s  =  1, 1  =  0,  and  j  =  1,  and 
final  state  is  3Pj,  i.e.,  s'  =  1,1'=  1,  and  summing  over  j'  =  0,  1,2, 
(C2(0I  I;00)  =  1  and  W2(01  lj';ll)  =  1/9),  we  get 

|e-  (n'W»>|2=T  £(2j'  +  l)|<nV'||i|nf)|2.  (A9) 

j' 

Substituting  in  equation  (A3),  we  get  the  equation  for  the 
radiative  lifetime 


P  =  T~ 


4  3 

m\c 


(3'£(2j'+i)i<nn'K>i2 


_1  (-Y  J- 

3 h  Vc/ 


^(2j'+l)|{nrHnf}|2 


£(2j'+D 


03|{n'fVK>|2. 


_4_  /w\3 
:  3h  Vc 


(A  10) 

giving  |(nT|n|n^)|  =  1.2  x  10~18  esu-cm.  Dasch  and 
Z?ec7ife/ [1981]  obtain  a  value  of  1.7  x  10“ 18  esu-cm  for  this 
matrix  element,  a  value  larger  by  factor  of  1.4  than  that 
calculated  here. 
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